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Abstract 



We show that density fluctuations in standard inflationary scenarios may take the 

pv^ . most general non-Gaussian distribution if the wave function of the Universe is not in 

^ \ the ground state. We adopt the Schrodinger picture to find a remarkable similarity 

^— N ' between the most general inflaton wavefunction and the Edgeworth expansion used 

cr^ . in probability theory. Hence we arrive at an explicit relation between the cumulants 

^^ \ of the density fluctuations and the amplitudes or occupation numbers of the various 

Qs ' energy eigenstates. For incoherent superpositions only even cumulants may be non- 

0^ ■ zero, but coherent superpositions may generate non-zero odd cumulants as well. 

(~| , Within this framework measurements of cumulants in Galaxy surveys directly map 

Qh! the wavefunction of the Universe. 

o 



1 Introduction 



Gaussian distributed density fluctuations are seen as a hallmark of the infla- 
tionary paradigm [1-4] . In addition, Gaussianity greatly simplifles calculations 
as only one set of parameters, the power spectrum, needs to be calculated. For 
these reasons, Gaussianity is commonly assumed as a matter of course. How- 
ever in non-minimal models of inflation [5-7] it is possible to generate non- 
Gaussian statistics. The magnitude [8] and type of departures from Gaussian- 
ity in these models appear to be restricted; in the simplest cases fluctuations 
are the square, or a power, of a Gaussian random fleld [7]. Furthermore, there 
have been a number of recent claims of detection of a non-Gaussian signal 
in the COBE data [9-15]. It is therefore of interest to explore further non- 
Gaussian inflationary models. For instance, one may wonder whether inflation 
can ever produce fluctuations with the most general probability distribution. 
It could be the case that some types of non Gaussianity are barred from even 
the most contrived inflationary models. 
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In this paper we consider standard inflationary models, and revisit the deriva- 
tion of the Gaussianity of density perturbations in these models. Microphysical 
fluctuations in the inflaton fleld then satisfy very approximately a linear har- 
monic oscillator equation. These should be subject to quantisation, with the 
common assumption being that the fleld lies in the ground state, which for an 
harmonic oscillator has a Gaussian wave function (eg. [16]). These microphys- 
ical fluctuations are then stretched by inflationary expansion, reaching what 
we nowadays call cosmological length scales. The square of their wavefunction 
provides their probability distribution, and since the square of a Gaussian is 
a Gaussian, these fluctuations have Gaussian statistics^. 

Let us now assume that inflation is triggered by quantum cosmology. However 
let us assume nothing from quantum cosmology (cf. [17,18]) about the initial 
form of the inflaton's wavefunction, ip, allowing it to be given by the most 
general superposition of energy eigenstates (a situation already considered in 
[19,20]). The wavefunctions ipn of these eigenstates are generally a Gaussian 
multiplied by a Hermite polynomial iJ„ [16]. The square of the wave func- 
tion, |?/^p, is the probability distribution of fluctuations emerging from such 
a scenario. The distribution we have just found has a striking resemblance 
to the Edgeworth expansion used in probability theory [21]. The Edgeworth 
expansion parameterises the most general distribution, and takes the form of 
a Gaussian multiplied by a series of Hermite polynomials. The coefficients in 
this series are polynomials in the cumulants of the distribution. Hence we may 
write a simple expression relating the components of ip and the cumulants of 
the density fluctuations' distribution. The purpose of this Letter is to estabhsh 
this interesting connection between observational and quantum cosmology. 

The plan of this Letter is as follows. In section 2 we review inflationary fluc- 
tuations, purposefully adopting the Schrodinger picture. In section 3 we re- 
view the other tool used in this Letter - the Edgeworth expansion. Then in 
sections 4 and 5 we establish the connection between the amplitudes of the 
various modes in the wavefunction and the cumulants of the distribution of 
the ensuing density fluctuations. We flnd that for an incoherent superposition 
of eigenmodes one may only obtain non-Gaussian distributions with non-zero 
even cumulants. However for coherent superpositions it is possible to generate 
the most general non-Gaussian distribution. Our results therefore show that if 
we take standard inflation for granted, we can in principle experimentally map 
the wavefunction of the Universe. The implications of such a measurement are 
briefly discussed in a closing section. 



^ Contrary to a popular myth, the Gaussianity of standard inflation's fluctuations 
is not due to the central limit theorem. 



2 Inflationary fluctuations in the Schrodinger picture 



We first review the formahsm used for describing inflationary fluctuations in 
the Schrodinger picture (eg. [19,22]). We assume that inflation arises from a 
single scalar field, the inflaton field 0. We write as an homogeneous back- 
ground plus a perturbation 5(j): 

0(x,t) = 0(t) + 50(x,t) (1) 



with 50 ^ 0. It is standard to decompose the perturbation 50(x, t) into 
Fourier modes 0k. In terms of these, the energy density fluctuations are given 

by [19,22]: 

(5p(k,t)=0(t)0k(t) (2) 



Density fluctuations are therefore proportional to 0k, which in turn is pro- 
portional to 0k. Therefore deviations from Gaussianity in 5p may be studied 
directly in terms of the statistics of 0k (though this is not true in the models 
considered in [23]). 

At early times we can adopt the free fleld theory limit and the modes 0k 
are then completely decoupled and behave as individual harmonic oscillators 
with a time- dependent mass m{t) [19,22]. The form of m{t) is dictated by the 
form of the potential V^(0) driving the inflaton fleld. In the slow-roll limit, and 
for the relevant wavenumbers, the mass m{t) is very approximately constant. 
We can therefore use the standard quantum mechanical treatment for the 
harmonic oscillator [16] to flnd the wave function for 0k (hereafter abbreviated 
as ^ = ^(0k))- 

In general -0 may be written as a superposition of energy eigenmodes: 

^ = H ^^^n (3) 

where n labels the energy spectrum E = huj{n + 1/2). The ipn take the form 



^„(0) = a^J--^je ^ (4) 



with normalisation flxing C„ as. 



" (2"n!v^(To)i/2 • ^^^ 



Here o"g is the variance associated with the (Gaussian) probabihty distribution 
for the ground state iV'oP- In other words (Tq is the power spectrum as derived 
in standard calculations in which the field is assumed to be in the ground 
state. (To is therefore related to the potential V{(f)) and to m according to the 
standard formulae [19,22]. We shall work with Hermite polynomials Hn{x) 
defined as 

H^{x) = {-ire^'£^e-^' (6) 



and normalised as 

J e-^"Hr,{x)Hm{x)dx = 2"7r'/%\6nm. (7) 



The most general probability density for the fiuctuations in is thus P = 
|-?/'P, with ip given by (3). The ground state (or "zero-point") fiuctuations 
are Gaussian, but any admixture with other states will be refiected in a non- 
Gaussian distribution function [19,20]. 

It is worth explaining the apparent discrepancy between our expressions and 
those in [20]. In [20] the authors regard each mode 0^ as a complex variable, 
and find a wave function for its modulus and phase, such that the probability 
distribution is most naturally expanded in Laguerre polynomials. In contrast, 
we consider the real and imaginary parts of (pk as independent real variables, 
and find separate wavefunctions for each. In this way we find an expansion in 
Hermite polynomials, allowing for a simpler connection with the Edgeworth 
expansion. The two treatments are equivalent and the probability densities 
obtained in [20] can be re-expressed in terms of Hermite polynomials, and vice 
versa. The results in this Letter still apply, although the algebra is significantly 
more tedious. 



3 The Edgeworth expansion 



We now review the Edgeworth expansion [21]. This expansion is a particular 
form of a class of series known as the Gram-Charlier Type A series. These 
are used to construct a distribution P{x) by means of its n-th order moments 
fin or cumulants /t„. These series are convergent under general conditions for 
P(x) [21,24]. 

For a standardized variable x (ie (x) = and (x^) = 1) the Edgeworth expan- 



sion takes the form 



-x^jl 



P(x) 



^CnHn{x) . 



If X is not standartized this expression can be easily modified. The coefficients 
Cn are tabulated polynomials in the moments /Xj or in the cumulants Kj. Un- 
fortunately the large body of complex formulae available in the literature uses 
Hermite polynomials Hn defined as 



Hn{x) = (-l)V 



n x^ /2_ 



dx^' 



-z2/2 



(9) 



as opposed to Hn{x) defined in (6). The two conventions may be bridged by 
means of [25] 



Hn{x/V2) 






(10) 



If P{x) is expanded around a Gaussian with the correct variance cr^ = (x^) 
(with (x) = 0) then the Edgeworth expansion is: 
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It is also possible to expand P{x) around a Gaussian with a different variance 
and mean. The coefficients in the expansion are then more complicated, but 
are listed in [21]. 



4 Coherent States 



Let us initially consider a coherent superposition of form (3). We first assume 
mild non-Gaussianity which we define through the condition |ao|^ 3> |ai|^, for 
i > 1. We can then consider the expansion of {ipl'^ only to first order in ctj, 



P(0) = 1^1^ 



e '"o 



27ro"n 



1 + E 



2^{ar. 



Hr, 



„Vj2"n!)V2 ny^^^ 



(12) 



where we have taken ao to have zero phase (so that to first order cto = !)• 

This leads to a one-to-one correspondence between the amphtudes of the var- 
ious energy eigenstates, and the combinations of cumulants appearing as co- 
efficients in the Edgeworth expansion. The latter simplify enormously if we 
only keep first order terms: we find that /t„ oc K(a„). Hence to first order 
the coherent contamination of the ground state by the n*^energy eigenstate 
is signalled by a non- vanishing cumulant /t„- For instance the presence of the 
third energy eigenstate results in n^ oc ^{a^) ^ 0, and, to first order, zero 
higher order cumulants. 

Of course we do not need to assume that |q;oP ^ Ictjp for i > 1, but then 
the series' convergence will have to be more carefully watched. For complete- 
ness, and for later use, we consider the more general case as well, making no 
assumptions about the relative amplitudes of the a^'s. We then have that. 



We may recover the Edgeworth expansion by noting that 



e-"'/J,(a;)/J,(a;)=e-"' 



(14) 



with. 



w = r _;L (IK) 

" {s - n)\{s - ty.{s - j)\ ^ ' 



with 2s = n + i+ j. This expression is subject to two selection rules: n + i + j 
must be even, and the triangular inequality \i— j\<n<i + j. One may 
derive (15) using (7) and the standard result for the integral over a product of 
three Hermite polynomials (formula 7.375.2 of [26]). Thus we obtain the more 
complicated expression 



^'■^^--^^'^'--vi^M^jm^yw <«) 



Hence the coherent admixture of the i*'^and j^'^eigenstates leads to terms in the 
Edgeworth expansion to all orders between \i—j\ and i+j. Note that in general 
one has (0) ^ 0, and that the variance of may be different from o"o, requiring 
suitable modifications to the coefficients in the Edgeworth expansion. 



5 Incoherent States 



If the various energy eigenstates are added incoherently we obtain the proba- 
bihty distribution, 



Pi<P) = I^P = E Kl'^Pl = ^ E !^^n( ^ ) (17) 

„ V2vrcTo „ 2"n! \V2aoJ 

As in the previous section this may be rewritten as 






P{<P) 



and one may now use the Edgeworth expansion results to infer the cumulants 
of the distribution. We see that in this case (0) = 0. However the variance of 
the perturbed distribution is never the same as o"q and is given by, 

a^ = alY,\anW2n+l) (19) 



a result which can be read off from (18) (with the aid of [21]), or easily derived 
in the Heisenberg picture. Hence formula (19), providing a^ — Uq, must be used 
when solving for the cumulants of the distribution from the coefficients in (18). 

The series (18) only contains Hermite polynomials of even order, since 6^" 
can only be non-zero for even k. Hence, we note the important result that 
incoherent superpositions generate non-Gaussian distributions with zero odd 
cumulants, in contrast with what was obtained in the previous section (as also 
mentioned in [19]). Furthermore for incoherent superpositions the presence of a 
non- vanishing n*^cumulant signals the presence of the n/2*'*energy eigenstate, 
in contrast with the result for coherent superpositions. In general, even to 
lowest order (which is second order in the amplitudes) one may not establish 
a one-to-one relation between cumulants /t„ and the l^np. The incoherent 
contamination of the ground state by the n'^'energy eigenstate generally leads 
to non-vanishing even cumulants of all orders up to 2n. 



6 Discussion 



We have shown that standard inflation may generate density fluctuations with 
the most general one-point distribution function, if we allow the inflaton fleld 



not to be in the ground state. According to this scenario the cumulants of the 
distribution provide a measurement of the components of the wave function 
of the inflaton field, in terms of energy eigenstates. The exphcit expression de- 
pends on whether the superposition is coherent or incoherent. In the coherent 
case all cumulants can be potentially non-zero. For incoherent superpositions 
only non-zero even cumulants can arise. 

There is however one type of non-Gaussianity which cannot be generated in 
this scenario: non-Gaussian inter-mode correlations. Translational invariance 
requires that the two-point function for modes 0k be diagonal ((0k0k') = 
Pfj,{k)6{k — 'k')). If the 0k are Gaussian they are fully determined by their two- 
point function, and so translational invariance requires them to be independent 
random variables. However if the reduced higher order correlators did indeed 
exist they would need not be diagonal as a result of translational invariance. 
For instance (0k</'k</'k) ^^^y t>e non-vanishing for all modes satisfying k = 
k' + k". Hence non-Gaussianity allows inter- mode correlations, even if we 
impose translational invariance. These correlation cannot be generated in the 
scenario we have proposed, because the modes 0k are genuine independent 
random variables. 

One may wonder whether the infiaton states we have considered are consistent 
with infiation. In most of the literature on infiationary density fiuctuations the 
inflaton is assumed to be in the ground state; two notable exceptions being 
[19] and [20]. This simplification has been justified using, amongst others, 
the rather dubious rationale that the field must be maximally symmetric, 
i.e. no scale should be privileged. However scenarios have been considered 
where a preferred scale is present, in particular that the Planck scale itself is a 
preferred scale [19]. A more physical argument for taking the ground state has 
also been proposed: that there must be a cutoff in the energy spectrum of the 
resulting inflaton particles to avoid inflnite energy densities [27]. The cutoff 
must also be low enough so that the energy density of the inflaton particles 
does not dominate over the potential driving inflation. The value of this cutoff 
is uncomfortably low if the inflaton particles were ever in thermal equilibrium. 
However it could happen that inflation is not preceded by a thermal phase, 
but starts straight out of quantum cosmology. In such a scenario one should 
not rule out mild admixtures of states above the ground state, with "mild" 
being deflned in the sense that any backreaction must be negligible. 

It is of course debatable whether or not such superpositions of states can 
be justifled by quantum cosmology. The Hartle and Hawking no-boundary 
proposal [17,18] and the Vilenkin tunneling wave proposal [28] imply that the 
inflaton wave function must be in the ground state. It would be interesting 
to investigate the implications to quantum cosmology of a detection of non- 
Gaussianity. 
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